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Abstract 

We discuss three specific modes of accretion disks 
around rotating magnetized neutron stars which may 
explain the separations of the kilo Hertz quasi peri- 
odic oscillations (QPO) seen in low mass X-ray bina- 
ries. The existence of these modes requires that there 
be a maximum in the angular velocity of the accreting 
material, and that the fluid is in stable, nearly cir- 
cular motion near this maximum rather than moving 
rapidly towards the star or out of the disk plane into 
funnel flows. It is presently not known if these condi- 
tions occur, but we are exploring this with 3D mag- 
nctohydrodynamic simulations and will report the re- 
sults elsewhere. The first mode is a corotation mode 
which is radially trapped in the vicinity of the max- 
imum of the disk rotation rate and is unstable. The 
second mode, relevant to relatively slowly rotating 
stars, is a magnetically driven eccentric (m = 1) os- 
cillation of the disk excited at a Lindblad radius in 
the vicinity of the maximum of the disk rotation. The 
third mode, relevant to rapidly rotating stars, is a 
magnetically coupled eccentric (m =1) and an ax- 
isymmctric (to = 0) radial disk perturbation which 
has an inner Lindblad radius also in the vicinity of 
the maximum of the disk rotation. We suggest that 
the first mode is associated with the upper QPO fre- 



quency, u u , the second with the lower QPO frequency, 
v& = v u — z/*, and the third with the lower QPO fre- 
quency, vg = u u — v tt /2, where is the star's rotation 
rate. 
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1 Introduction 

Low mass X-ray binaries often display twin kilo- 
Hertz quasi-periodic oscillations (QPOs) in their X- 
ray emissions (van der Klis 2006; Zhang et al. 2006). 
A wide variety of different models have been pro- 
posed to explain the origin and correlations of the 
different QPOs. These include the beat frequency 
model (Miller, Lamb, & Psaltis 1998; Lamb & Miller 
2001; Lamb & Miller 2003), the relativistic precession 
model (Stella & Vietri 1999), the Alfven wave model 
(Zhang 2004), and warped disk models (Shirakawa & 
Lai 2002; Kato 2004). 

A puzzling aspect of the some of the twin QPO 
sources considered in this work is that the differ- 
ence between the upper v u and lower v& QPO fre- 
quencies is roughly either the spin frequency of the 
star i/» (3 cases where = 270, 330, & 363 Hz) 
or one-half this frequency, v*/2 (4 cases where = 
401, 524, 581, & 619 Hz), for the cases where v* 



is known, even though v u and vg vary significantly 
(see, e.g., Zhang et al. 2006). A further type of be- 
havior appears in the source Cir X-l (Boutloukos et 
al. 2006), but this is not considered here. The cases 
where v u — « may be explained by the beat 
frequency model (Miller et al. 1998), but the expla- 
nation of the cases where v u — vi w v*/2 is obscure. 

Section 2.1 discusses the corotation instability, §2.2 
the execentric (m = 1) mode of the disk driven by the 
star's rotating magnetic field, and §2.3 the coupled 
execentric plus axisymmetric mode (m = & 1) also 
due to the star's rotating magnetic field. Section 4 
gives the conclusions. 

2 Three Modes 

2.1 Corotation Instability 

We assume a pseudo-Newtonian gravitational poten- 
tial $ g — -GM t /(r — rs), where M* is the star's 
mass and rs = 2GM»/c 2 . In the absence of the 
star the angular velocity of disk matter is Q g ^, — 
{GM*/[r(r-r s ) 2 ]} 1/2 for r > 3r s . Near the star the 
disk's angular rotation rate in the equatorial plane is 
modeled as 



1 + f(r) 1 + f(r) ' 



(1) 



where f(r) — exp[— (r — r )/A] with r the stand- 
off distance of the boundary layer and A its thick- 
ness which are expected to depend on the accre- 
tion rate and the star's magnetic field. The ra- 
dial force equilibrium in the midplane of an axisym- 
metric disk is /w(f2 2 - ) _ ^(p + B 2 /8n)/dr, 
where p is the midplane density and midplane field 
B = B(r)z. Figure 1 shows the equilibrium quan- 
tities for an illustrative case. Clearly, fl^,(r) has a 
maximum value outside of the star at a distance 
denoted r m ~ r . The importance of this maxi- 
mum for models of QPOs was discussed earlier by 
Alpar and Psaltis (2005). Three dimensional mag- 
netohydrodynamic (MHD) simulations of disk accre- 
tion to rotating magnetized stars can in principle 
be used to determine f2^(r) for different conditions 
(Romanova, Kulkarni, & Lovelace 2007). However, 



for the present purposes the dependence of equa- 
tion (1) is used. It is needed only for distances 
r > r m as discussed below. In this region the radial 

epicyclic frequency is f2 r (r) = [r~ 3 d(r 4 0, 2 i )/dr] 1 / 2 . 
The maximum value of O^(r) has the approximate 
dependence max(^/ 27r ) ~ 2040(3r s /r ) 1 ' 74 Hz for 
3 < r /r s < 5, A/r s = 0.05, and M* = 1.4M . For 
this range of r , max(f2^,) changes by a factor of 2.4. 
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Figure 1: Equilibrium profiles of the main variables 
normalized to their maximum values. For this case, 
M* = 1.4M , = VL^/2tt = 300 Hz, r s w 4.14 x 10 5 
cm, r = 4r s , A = 0.05r s , c s = O.lrO^, p max ~ 
0.024 g/cm 3 , B(r/r s = 3.5) 
max(fJ /27r) w 1220 Hz. 



1.7 x 10 9 G, and 



We consider a WKB treatment of the corotation 
(u> pa mfi^) or Rossby type wave of the disk with 
pressure perturbation 



bp ~ exp 



drk + im<f) — iiot 



where k is the radial wavenumber, m = 1,2, .., and 
u> = u> r + iu>i, with u r the angular frequency of the 
perturbation and u>i the growth rate. For the condi- 
tions of Figure 1, the Alfven speed ca — B/(4iTp) 1 ^ 2 
is much larger than the sound speed c s - Also, we as- 
sume and verify later that \u> — mVt<p\ 2 <C ft 2 . Under 
these conditions 
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where Alj = cj — mfioi, L] 



Alu 2 
(2) 

1 = dln(B/p)/dr, L- 1 = 



d]n(B)/dr, T = pQ^/Q 2 , which all depend on r, and 
SR(..) denotes the real part (Lovelace, Turner, & Ro- 
manova 2007). Equation (2) generalizes the calcula- 
tion of Lovelace et al. (1999) to include the influence 
of the magnetic field perturbation <5B = 5B(r, <f>,t)z. 

Figure 2 shows the radial dependence of k 2 for a 
representative case. For the chosen value w r /27r = 
1100 Hz, which is somewhat less than the maximum 
value of r^(r), and m = 1, it is seen that k 2 (r) > in 
a finite radial interval in the vicinity of the maximum 
of Q$(r). Thus the wave is radially trapped in the 
vicinity of the maximum of fi^r). Analogous radi- 
ally trapped modes were analyzed earlier by Lovelace 
et al. (1999) and Li et al. (2000) and verified in 
two dimensional hydrodynamic simulations (Li et al. 
2001). The Bohr-Sommerfeld quantization condition 



Jr(in) ^ r k = + 1/2)tt, n = 0, 1, .. allows the deter- 
mination of the growth rate uii. For the case shown, 
U!i/2TT — Vi w 55 Hz for n = which gives the largest 
growth rate, and r(m)/rs = 4.08, r (out) /rs = 4.48. 
The growth rate increases as A decreases. Similar 
values of the growth rates are found for u r somewhat 
less than mmax(f20) for m = 2, 3, ... The nonlinear 
saturation of the growth of the modes can in prin- 
ciple be found by MHD simulations (Koldoba et al. 
2002; Romanova et al. 2007). 

From Figure 2 we see that the validity of equa- 
tion (1) is needed only from the vicinity of the max- 
imum of fi^(r), that is, from rjr$ — 4.08 where 
fi^/max(fi^) = 0.94 (the inner turning point) and 
outward (including the outer turning point at r /r$ — 
4.48). 

Owing to the perturbation, the surface tempera- 
ture of the disk is T(r, <p, t) — 



To + 6Ti exp(i0 — u>it) +ST2(r) exp(2i0 — iu>2t) + ■ 



where T)(r) is the unperturbed temperature, 
STi : 2(r) <C To are the amplitudes of the m = 1,2 
corotation modes, and are their frequencies. 

The corresponding flux density is proportional to 
S(r,<f>,t) ~ T 4 + 4T 3 JT 1 exp(i</> - wrf) + .. The to- 
al flux for a face-on disk, L ~ J rdrdipS, is indepen- 
ent of time. For a more general disk orientation, the 
Doppler effect due to the disk rotation gives a boost 
for say 4> — and a decrement for cj> = vr. This cor- 
responds to multiplying S by [1 + eexp(— i<f>)], with 
e(r) -C 1. Consequently, there is a contribution to 
the source flux SL ~ J rdrdtfie~x.p(— i<fi)S(r, (ft, t) ~ 
4 J rdrT^ eSTi(r) exp(— iuj\t). As explained in the 
next section, we interpret this frequency as the upper 
frequency component of the twin QPOs. 
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Figure 2: Radial dependences of the different fre- 
quencies Vfy — Q$/2ir, the resonant frequency v res = 
u} r /2w — 1100 Hz, and the square of the radial 
wavenumber, k 2 obtained from equation (2) for m = 
1 and for the same conditions as for Figure 1. 
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2.2 Magnetically Driven m — 1 Mode 

We now give a heuristic treatment of perturba- 
tions of the disk excited by the star's rotating, non- 
axisymmetric magnetic field. One component of a 
general perturbation is described by the radial dis- 
placement of the disk matter, £(r,<f>,t) (e.g., Zhang 
& Lovelace 2005) The equation for £ is simply 

^ = -(^ + ^)e + 6F(r,4>,t), (4) 

where d/dt = d/dt + fl^d/drf*, and SF is the radial 
force due to the star's magnetic field. The radial os- 
cillation frequency (squared) on the right-hand side of 
the equation consists £l 2 as given earlier plus a ther- 
mal contribution k 2 c 2 , where k is the radial wavenum- 
ber of the perturbation in the WKB approximation. 

We consider the case where the star has a small 
size or point-like dipole located close to the star's 
surface near the rotation axis with its magnetic mo- 
ment parallel to the star's surface. Ruderman (2006) 
gives physical arguments for this type of configura- 
tion of spun-up neutron stars. At a given radial dis- 
tance in the equatorial plane, the total magnetic field 
B z — B z + B\ consists of the vacuum component 
B z = B (r) cos(0 — f2*i), with fi» the star's angular 
rotation rate, and the induced component, B\. The 
component B v z acts to drive a current flow in the 
disk /C<£ which in turn gives rise to the induced mag- 
netic field of the disk B\. The radial force includes 
a coherent component cx B v z which can be written as 
5F = dexp(i<f)-in t t), where d{r) = (/C^)S /(cS) 
with the average is over <p. Considering only the co- 
herent component, equation (4) is then a driven os- 
cillator. With ik — > d/dr and £ <~ exp(i(f> — iQ*t), it 
becomes 

(^- p w)^ = - C i/ c '' ( 5 ) 

with V = [Q 2 — (f^ — VL^) 2 ]/c 2 . Depending mainly 
on the value of Q*, there may be Lindblad reso- 
nances with inner and outer Lindblad radii r Li ^ a 
where T>(rLi,o) — 0. The region between the two 
radii is permitted in the sense that T> < while the 
regions r > and r < ru are forbidden. The exis- 
tence of the Lindblad resonances means that a weak 



magnetic disturbance can give rise to a strong disk 
response in the vicinity of tl proportional to Ci[tl) 
(Goldreich & Tremaine 1979). Because of the rapid 
decrease of the magnetic field the response at the 
inner Lindblad resonance is expected to be stronger 
than that at the outer resonance. However, the solu- 
tion of equation (5) is beyond the scope of this work. 
We find that there are Lindblad resonances only for 
for z/* = £!*/27r < 380 Hz for the same conditions as 
Figure 2. For higher z/„, all values of r are forbidden. 
For = 300 Hz, r^i/rs ~ 4.34, which is within the 
region of the mentioned corotation instability, and 
fLo/fS ~ 4.83 which is outside the region of corota- 
tion instability. The disk velocity fi^r is supersonic 
relative to the velocity of the perturbation f2*r for r 
larger than rn- Consequently the perturbation is a 
leading spiral wave with k < 0. At the outer Lindblad 
radius the disk velocity is also supersonic relative to 
0*r, and excitation of the disk at this radius gives a 
trailing spiral, k > for r < r^ Q . 

The interaction between the corotation pertur- 
bation and the magnetic perturbation of the disk 
is in general nonlinear. The perturbed disk sur- 
face temperature can be represented as a prod- 
uct of the two perturbations T(r, <fi, t) ~ [1 + 
€M\£/r\ exp(— i(f> + iQ*t)][To + <57\ exp(i</> — iw\t) + 
..} = T + e M \£/r\6Ti exp[-i(wi - + .., where 
EM < 1. Consequently, there is a contribution to 
the source flux SL <~ f rdrd<f>exp(—i(j))S(r,<j>,t) ~ 
4 / rdrT 3 e M |£/r|<5Ti(r) exp[-i(wi-fi*)t], where S ~ 
T 4 . For <^ 380 Hz, we interpret u>\ — fi* as the 
lower frequency component of twin QPOs. 

2.3 Magnetically Coupled m = k 1 
Mode 

Consider now higher rotation frequencies v*. Note 
that the radial force perturbation includes a con- 
tribution of the form &F ~ Yr Y %t(£ )d{JC 4> B z / c) / dr 
which again has a coherent term proportional to 
B z . We consider only the coherent term which can 
be written as 6F = $t(£ )D\ cos(</> — Q*i), where 
Di = d{JC<pB ) / dr / (cY.) — real. In this case £ neces- 
sarily consists of different 0— harmonics. That is, £ = 
£o exp(—iu;t)+£i cxp(i<f> — iut) + £2 exp(2i<fi — u)t) + ... 
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We find 

(f s ^i-W-{^-u,f] s je 1 = -^1(6) 

where necessarily u> = f2»/2. For the radii of interest 
in the vicinity of r , ft 2 , ^> uj 2 , so that the right- 
hand side of the second equation is w — -Df£i/(4f2^). 
We identify |-Di(r)| as the Alfven frequency squared 
Q\ at the distance r and estimate that Q 2 A <C f^- 
Evidently the equation for £ 1 corresponds to free os- 
cillations. It has Lindblad radii approximately where 
V{r) = Q 2 r - (f^ - uj) 2 = 0. For u* = lu/tt = 600 
Hz and the other parameters the same as in Fig- 
ure 2, the inner Lindblad radius is at r^j/rs = 4.2 
while the outer Lindblad radius is at a large distance 
(r/rg > 6). This mode may be driven at rn by 
noise or fluctuations in the disk at this radius, and 
it is a leading spiral wave. The mode amplitude can 
in principle be found using three-dimensional MHD 
simulations (Koldoba, et al. 2002; Romanova et al. 
2007). 

The influence of the magnetically coupled modes 
on the flux follows the discussion of §2.3. 
We find SL ~ frdrd(f>exp(—i<f))S(r,(j),t) ~ 
4 / rdrT§ e M \£i/r\6Ti(r) exp [-£ - Q*/2)t]. For 
>^ 380 Hz, we interpret u)% — f2»/2 as the lower 
frequency component of twin QPOs. 

3 Conclusions 

We discuss three modes of accretion disks around 
rotating magnetized neutron stars which may ex- 
plain the frequency separations of the twin kilo-Hertz 
QPOs seen in accreting X-ray binaries. The existence 
of these modes requires that there be a maximum 
in the angular velocity of the accreting material and 
that the fluid is in stable, nearly circular motion near 
this maximum rather than moving rapidly towards 
the star or out of the disk plane into funnel flows. 
It is presently not known if these conditions occur, 
but we are exploring this with 3D magnetohydrody- 
namic simulations and will report the results else- 
where. The first mode is a corotation mode which 



is radially trapped in the vicinity of the maximum 
of the disk rotation rate and is unstable. A sim- 
ple dependence is assumed for the angular rotation 
rate of the disk fta (r) which has a maximum at a ra- 
dius r m outside the neutron star. The unstable mode 
has a frequency lo\ somewhat less than max^^,) and 
this can vary by a significant factor depending on 
the state of the disk (e.g., the accretion rate). We 
suggest that this mode is associated with the upper 
QPO frequency, v u = lui/2tt. The second mode is a 
magnetically driven eccentric (m = 1) oscillation of 
the disk excited at the inner Lindblad radius which 
is in the vicinity of the maximum of the disk rota- 
tion. The star's magnetic field is assumed to be the 
form discussed by Ruderman (2006). The Lindblad 
radii occur only for relatively slowly rotating stars, 
<^ 380 Hz. For these stars we suggest that the 
lower QPO frequency is vg = v u — v*. The third 
mode, relevant to more rapidly rotating stars, is a 
magnetically coupled eccentric (m = 1) and an ax- 
isymmetric (m = 0) radial disk perturbation. It has 
an inner Lindblad radius also in the vicinity of the 
maximum of the disk rotation. For these stars the 
lower QPO frequency is Vu = v u — v*/2. A problem 
remaining for future work is the determination of the 
saturation amplitudes of the different modes. 
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